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Constraints on dark energy from H II starburst galaxy apparent 

magnitude versus redshift data 

Data Mania^'^ and Bharat Ratra^ 
ABSTRACT 



In this pap er we use H I I star burst galaxy apparent magnitude versus red- 
shift data from ISiegel et al.l ( 120051 ) to constrain dark energy cosmological model 
parameters. These constraints are generally consistent with those derived us- 
ing other data sets, but are not as restrictive as the tightest currently available 
constraints. 



1. Introduction 



There is significant observational evidence supporting the idea that we live in a spatially- 
fiat Universe that is currently undergoing accelerated cosmological expansion. Most cosmol- 
ogists believe that this accelerated expansion is powered by dark ene rgy which dominates the 
current cosmological energy budget (for reviews of dark energy see iTsujikawal l2010l : iLinder 
2010l : iBlanchardI l2010l : ISapond l2010l . and references therein) o 



T he energy bud get of the "standard" model of cosmology — the spatially-flat ACDM 
model (jPeeblesI 119841 ) — is currently dominated by far by dark energy: Einstein's cosmolog- 
ical constant A contributes more than 70 % of the budget. Nonrelativistic cold dark matter 
(CDM) is the next largest contributor (more than 20 %), followed by nonrelativistic baryons 



^Department of Physics, Kansas State University, 116 Cardwell Hall, Manhattan, KS 66506, USA mania, 
ratra@phys.ksu.edu 



gia 



^Center for Elementary Particle Physics, Ilia State University, 3-5 Cholokashvili Ave., Tbilisi 0179, Geor- 

i 

^ Some, however, feel that this accelerated expansion is better v iewed as an indication that the general 



relatiyistic description of g ravitation needs to be improved on (see Ide Felice fc Tsuiikawa 



2010; 



Tsuiikawa 



2OIOI : IJain fc Khourvll2010l and references therein) . In this paper we assume that general relativity provides 
an adequate description of gravitation on cosmological scales. 



-2- 



(around 5 %). For a review of the standard model see lRatra fc Vogeleyl (120081 ) and references 
therein. It has been know for some time t hat the ACDM model is reasonably consistent with 
most observational constraints 



see, e.g. 



200?! : kUen et al.l 120081 for early indications)^ 



Jassal et al. 2010: Wilson et al. 2006: Davis et al. 



However, the ACDM model appears to leave some conceptual questions unanswered. 
One of these is that the naive energy density scale we expect from quantum field theory 
considerations is many orders of magnitude higher than the measured cosmological constant 
energy density scale. Another is what is known as the coincidence puzzle: The energy density 
of a cosmological constant does not change over time but the matter density decreases with 
the cosmological expansion, so it is puzzle why we observers live at this seemingly special 
time, when the dark energy and nonrelativistic matter densities happen to be of comparable 
size. 

These puzzles, and perhaps others, could be alleviated if we assume that the dark energy 
density was higher in the past and slowly decrease d in time, thus remain ing comparable to 



the nonrelativistic matter density for a longer tim ejRatra fc PeebleslEoSsI ) . Many such time- 
varying dark energy models have been proposedjj In this paper, for illustrative purposes, 
we consider two dark energy models and one dark energy parametrization. 

In the ACDM model, time-independent dark energy — the cosmological constant — 
may be modeled as a spatially homogeneous fluid with equation of state parameter cja = 
Pa/Pa = —1 (where and pA are the fluid pressure and energy density). 

In an attempt to describe slowly decreasing dark energy, the dark energy fluid equation 
of state is extended to the XCDM X-fiuid parametrization. Here, dark energy is again 
modeled as a spatially homogeneous (X) fluid, but now with an equation of state parameter 
wx = Vxl Pxi where wx[< —1/3) is an arbitrary constant and "px and px are the pressure 
and energy density of the X-fluid. When wx = — 1 the XCDM parametrization reduces 
to the complete and consistent ACDM model. For any other value of wx[^< —1/3), the 
XCDM parametrization is incompl ete as it cannot describe spatial inhomogeneities (see, e.g. 
Ratralll991t iPodariu fc Ratrall2000[ ) . For computational simplicity, here we study the XCDM 
parametrization in only the spatially-flat cosmological case. 



A slowly rolling scalar field 0, with decreasing (in 0) potential energy density V^(0), can 



The ACDM model assumes the "standard" CDM structure formation pict ure, which might be in some 



observational difficulty (see, e.g 
^ For recent discussions see, e. 



Peebles fc Ratralboost fPerivolaropouloslboid ). 



Pettorino et al.l (|2010l ). iLiu et al 



Noyosvadlvi et al (120101) .iHarko fc Lobol (120101 ) . iKeresztes et al.l (|2010f ). 
(|2010l) . ICostal (|2010l ) iFaraiollahi et al.l (|201l[ ) and references therein. 
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be used to consistently model a gradually decreasing, time-evolving, dark energy density. An 
inverse power-law potential energy density V{(j)) oc ( ^~", where a is a non- negative constant, 
provi des a simple realization of this (pCDM scenario ([Peebles fc Ratralll988l : iRatra fc Peebles 
19881 ). When a = the 0CDM model reduces to the corresponding ACDM case. For com- 
putational simplicity, here we study the 0CDM model in only the spatially-flat cosmological 



case. 



Current observational data convincingly indicate that the cosmological expansion is ac- 
celerating. The main support for accelerated expansion comes from three types of data: 



AmanuUah et al.ll2010t iHolsclaw et al.ll2010l : IWang et aL 



crowaye background 



2009: Liu et al. 2011 



Wane et al. 


2011; 


Guv et al. 


^^^^^^^^^ , v^.^., 

2010): cosmic mi- 


fsee, e.g., 


Podariu et al. 


l2001b: 


Komatsu et al. 



La Vacca fc Bonomettol 2011 ) combin ed with low estimates of the 



cosmological matter density (see, e.g, IChen fc Ratral l2003bl): and, baryon acoustic oscil 



lation (BAO) peak length scale estim ates (see, e.g., ISamushia fc Ratral l2009t IWangI 12010 



Percival et al.ll2010l : iBlake et al.ll201ll ). However, the error bars associated with these three 
kinds of data are still too large to allow for a significant observational discrimination between 
the ACDM model and the two simple time-varying dark energy models discussed above. 

There are two major reasons to consider additional data sets. First, it is of great 
interest to check the above results by comparing to results derived from other data. If the 
new constraints differ significantly from the old ones, this could mean that the model being 
tested is observationally inconsistent, or it could mean that one of the data sets had an 
undetected systematic error. Either of these is a significant result. On the other hand, if 
the constraints from the new and the old data are consistent, then a joint analysis of all the 
data could result in significantly tighter constraints, and so possibly result in significantly 
discriminating between constant and time-varying dark energy. 

Other data that have recently been us ed to constrain da r k energy raodels i nclude strong 



gravitational lensing mea surements (e.g., IChae et al.l l2004l : iLee Sz Na 120071 : IZhang &: Wu 



2010 



2000 



Biesiada et al 



Chen fc Ratra 



20101). angular size as a function of redshift observations (e.g.. lGuerra et al. 



shift measurements (e.g.. 



Samushia et al. 



2008 



2007: 



Pan et al. 



Samushia fc Ratra 



ture observations (e.g . 



Brouzakis et al.ll2010t ICampanelli et al. 



2008 



2003at iBonamente et al.l 20061). Hubble paramete r as a function of red- 



Jimenez et al. 20031: Samushia &: Ratra 2006: Sen &: Scherrer 



2010h. galaxy cluster gas mass fractio n data (e.g. 



Ettori et a. 



Baldi fc Pettorino 



20081: 



Allen et al. 



2009 



2011 



20111 



iTong fc Nohl I2OIII ) , and large-scale struc- 



De Boni et al. 2011 



Mortonson et al. 2011 



While constraints from these data are less restrictive than those derived from the SNela, 
CMB and BAO data, both types of data result in largely compatible parameter restrictions 
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that generally support a currently accelerating cosmological expansion. This gives us con- 
fidence that the broad outlines of the "standard" cosmological model are now in place. 
However, there is still ambiguity. Current observations are unable to differentiate between 
different dark energy models. For instance, while current data favor a time-independent 
cosmological constant, they are unable to rule out time- varying dark energy. More and 
higher-quality data is needed to accomplish this important task. 

It is anticipated that future space missions will result in significantly more and better 



SNela, BAQ, and CMB anisotropy d ata (see, e.g.. |Podariu et al.ll2001al : ISamushia et al.ll2011 



Wane et al. 


2010; 


Astier et al. 


2010) 



tests that make use of different sets of objects . Recent examples include the lookback time 
test (e.g., iPires et al.l 2006t iDantas et al.ll201l[) and the gamrna-ray burst luminosity versus 



Schaefer 



2007 



Wang||2008l : IXu fc Wanal2010l ). Gamma-ray bursts, in 



redshft test (see, e.g. 

particular, are very luminous and can be seen to much higher redshift than the SNela and 
so probe an earlier cosmological epoch. H II starburst galaxies are also very luminous and 
can be seen to redshifts exceeding 3. Recent work has indicated that H II starbusrt galaxie s 
might be standardizable candles (IMelnick et al.ll2000t iTerlevich et al.l 120021 : lMelnickll2003l ). 
because of the correlation between their velocity dispersion, H fj luminosity, and metallicity 
jMelnicklllQTsI : IXerlevich fc Melnickl EoSli iMelnick et aPllQSsh . 



In this paper we use H II galaxy data from lSiegel et al.l (2005 ) to constrain param eters 
of the th ree dark energy ni odels mentioned above. iPlionis et al.l ( 120091 . l2010l 1201 ll ) have 
used the ISiegel et al.l (|2005[ ) data to constrain the XCDM parametrization. Here we also 
constrain parameters of the consistent and physically motivated ACDM and 0CDM cosmo- 
logical models. We also derive constraints on the parameters of t hese models and the XCDM 
parametrization from a joint analysis of the ISiegel et al.l (120051 ) H II galaxy data and the 
Percival et al.l (120101 ) BAO peak length scale measurements. 



In the next section we list the basic equations for the dark energy models we consider. 
In Sec. [3] we discuss the analysis method used for the H II galaxy data and the resulting 
constraints and how they compare with those from other cosmological tests. Section H] 
presents constraints from a joint analysis of H II and BAO data. Section |5] summarizes our 
results. 



2. Dark energy model equations 

The Friedman equation for the ACDM model with spatial curvature is 
H{z, Ho, p) = Hoy/nUl + z)^ + I^a + 1)r(1 + z)^. 



(1) 
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Here Hq is the Hubble constant and H{z) is the Hubble parameter at redshift z, flm, ^a, 
and = 1 — Qm — are the nonrelativistic (baryonic and cold dark) matter, cosmolog- 
ical constant, and space curvature energy density parameters, respectively, and the model 
parameter set p = (fim, ^a)- 

For the spatially-flat XCDM parametrization, the Hubble parameter is 



H{Z, Ho, p) = Ho^n^il + Z)^ + (1 - + ;2)3{l+-x). (2) 

Here wx < —1/3 and the model parameter set p = {Qrn,<^x)- 

For the spatially-flat 0CDM case we have a coupled set of equations governing the scalar 
fleld and Hubble parameter evolution 

+ 3i/0 - ("^^^ = 0, (3) 



H{z, Ho, p) = HoJilmil + + n^{z). (4) 



Here we consider the inverse power law scalar fleld (0) potential energy density (jPeebles fc Ratra 



19881) so 



where G is the Newtonian gravitational constant and a > is a free parameter (that deter- 
mines k). In this case the model parameter set is p = {Qm, «)• 

The distance modulus function of redshift is 

Ho, p) = 25 + 5 log [3000(1 + z)y{z)] - 5 log h, (6) 

where h is the dimensionless Hubble constant deflned from Ho = lOOh km s~^ Mpc~^. In 
the spatially open case y{z) is 



y{z, P) = -j== sinh / T^TTT ' (7) 




where E{z) = H{z)/Ho. When Qr = this reduces to y{z) = J dz'/E{z'). The angular 



diameter distance cIa = cy{z)/[Ho{l + z)]. 
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3. Constraints from H II galaxy apparent magnitude data 



We use the 13 yUobs(-2j) measurements of (ISiegel et al.ll2005l ). listed in ITable T| to con- 
strain cosmological parameters by minimizing 



XhiiI^cPJ = -2 ■ 



Here Zi is the redshift at which /iobs(-2j) is measured, /ipred(-Zj, Hq, p) is the predicted distance 
modulus at the same redshift for the model under consideration, and cTj is the average of the 
upper and lower error bars listed in ITable II 



The ISiegel et al.l (120051 ) measurements listed in ITable II are computed from 



/iobs = 2.5 log 



a 



HI3 



2.5 log 



AHf3 + Zq 



(9) 



where FHFi and Ah h are the Hp flux and extinction and 0/H is a metallicity. Following 
Plionis et al.l ( 120111 ). for the zero point magnitude we use Zq = —26.60, use Hubble constant 
value Hq = 73 km s~^ Mpc~^ (and do not account for the associated uncertainty), and 
also exclude two H II galaxies (Q1700-MD 103 and SSA22 a-MD41) that show signs of a 
considerable rotational velocity component ( lErb et al.ll2006l ). 



^ijj. The la, 2a, and 3a 



The best flt parameter set p=K is where Xmi is at a minimum, x- 
confldence intervals, in the two-dimensional parameter space, are enclosed by the contours 
where = xLin + with Ax^ = 2.30, Ax^ = 6.17, and Ax^ = 11-8, respectively. 



The H II galaxy data constraint s on cosmolog i cal p arameters are shown in Figs. 1 — 
Our results generally agree with iPlionis et al.l ( 120111 ) for the XCDM parametrization 



(compare our Fig. 2 and their Fig. 10). The small differences are due to the fact that in 
our analysis we do not use specially weighted sigmas, but rather average distance moduli 
uncertainties, and also ignore gravitational lensing effects. Both analyses ignore small Hq 
uncertainties; these are not important for our illustrative purposes here, but should be 
accounted for in an analysis of improved near-future H II galaxy data. 

The H II galaxy data constraints in Figs. 1 — 3 are not as restrictive as those that 
follow from SNela, BAO, or CMB anistro py data. They are, ho wever, comparable to those 
from Hubble parameter observ ations (see IChen fc Ratrall2011bl . and references therein) or 
lookback time observations (see lSamushia et al.ll2010l . and referen ces therein) , and some what 
more restrictive than angular diameter distance constraints (see Chen fc Ratra boila, and 



references therein) and gamma-ray burst luminosity distance ones (see ISamushia fc Ratra 



2010l . and references therein). We again note that we have not accounted for the uncertainty 
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in Ho in our analysis, so the H II galaxy constraints derived here are a little more restrictive 
than they really are. However, our analysis clearly illustrates the potential constraining 
power of near-future H II galaxy luminosity distance data. 



4. Joint constraints from BAO and H II galaxy data 



To constrain observables using BAO dat a we follow the proc edure of iPercival et al. 



teoioh . With Dv{z) = [(1 + z)^dlcz/H{z)Y/^ 
£1^(0.275) = (1104 ±30) 



Percival et al.l (120101 ) measure 

-0.134 / ^ , 2 \ -0.255 



0.02273 



0.1326 



Mpc. 



(10) 



We use this to build the likehhood estimator £bao with a Gaussian prior of fl rnh"^ = 0.1326± 



0.006 3, and neglect the error for as WMAP5 data constrains it to 0.5 % (IKomatsu et al. 
20091 ). The BAO data constraint contours are shown in Figs. 4 — 6. 



To derive joint H II galaxy and BAO data constraints we maximize the products of 
likelihoods C{p) = £hiI'^^bao! where £hii = e^^Hn/^^ to get the best fit set of parameters p*, 
and the 1, 2, and 3a contours defined as points where the likelihood equals e~^-^*^/^, e"^'^''/^, 
and e"^^'*/^, respectively, of the maximum likelihood value. These contours are shown in 
Figs. 4 — 6. Clearly, even adding currently available H II galaxy data to the mix tightens the 
constraints. 



Conclusion 



We have used the starburst galaxy luminosity distance data of ISiegel et al.l (|2005[ ) to 
constrain cosmological parameters in some dark energy models. The resulting constraints 
are largely consistent with those that follow from other currently available data although 
they are not as restrictive as the constraints derived from SNela, BAO, and CMB anisotropy 
data. 



The ISiegel et al.l ( l2005l ) H II data are preliminary H II data. It is gratifying that they 
result in interesting constraints on cosmological parameters. We anticipate that new, near- 
future, H II galaxy data will provide very useful constraints on cosmological parameters that 
will complement those derived using other data. 



DM thanks Manolis Plionis, Glenn Horton-Smith, and Yun Chen for valuable sugges- 
tions. We acknowledge support from DOE grant DEFG030-99EP41093, NSF grant AST- 
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Fig. 1. — H II galaxy data 1, 2, and 3a confidence level contours in the (fi^, ^^a) plane for 
the ACDM model. The dotted line corresponds to the spatially-flat ACDM case and the 
shaded area in the upper left hand corner is the part of parameter space without a big bang. 
The best-fit point with Xmm — 53.3 is indicated by the solid black circle at Qm = 0.19 and 
= 0.98. 



'0.0 0.1 0.2 0.3 0.4 0.5 

Fig. 2. — H II galaxy data 1, 2, and 3(7 confidence level contours in the (O^, wx) plane for 
the spatially-flat XCDM parametrization. The dotted hne corresponds to the spatially-flat 
ACDM case. The best-fit point with Xmin — ^3.3 is indicated by the sohd black circle at 
Q„ = 0.17 and wx = -0.86. 
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Fig. 3. — H II galaxy data 1, 2, and 3a confidence level contours in the {flm, a) plane for the 
spatially-flat (f)CDM model, a = corresponds to the spatially-flat ACDM case. The best-fit 
point with xLin — ^3.3 is indicated by the solid black circle at fl^ — 0.17 and a = 0.39. 
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Fig. 4. — Joint H II galaxy and BAO data (solid lines) and BAO data only (dashed lines) 1, 
2, and 3a confidence level contours in the {flm, ^a) plane for the ACDM model. Conventions 
and notation are as in Fig. 1. The best-fit point with — 21og(£niax) = 55.2 is indicated by 
the solid black circle at Qm = 0.25 and Qa = 0.95. 
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0.15 0.20 0.25 0.30 0.35 0.40 

Fig. 5. — Joint H II galaxy and BAO data (solid lines) and BAO data only (dashed 
lines) 1, 2, and 3a confidence level contours in the {^m, ^x) plane for the spatially-flat 
XCDM parametrization. Conventions and notation are as in Fig. 2. The best-fit point with 
—2 log (Lmax) = 53.5 is indicated by the solid black circle at flm — 0.25 and Wx — —1.41. 
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Fig. 6. — Joint H II galaxy and BAO data (solid lines) and BAO data only (dashed lines) 1, 
2, and 3a confidence level contours in the {Qm, cx) plane for the spatially-flat 0CDM model 
Conventions and notation are as in Fig. 3. The best-fit point with —2 log(>Cinax) = 55.6 is 
indicated by the sohd black circle at Jl^ = 0.27 and a — 0. 
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Table 1: ISiegel et al.l ( 120051 ) H II starburst galaxy distance moduli and uncertainties 



Galaxy 


z 


Atobs ± Cr 


Q0201-B13 


2.17 


47.4911111 


Q1623-BX432 


2.18 


45.451^:^? 


Q1623-MD107 


2.54 


44.82l?:i 


Q1700-BX717 


2.44 


46.641?:^^ 


CDFa CI 


3.11 


45.77l?:i 


Q0347-383 C5 


3.23 


47.121°:^^ 


B2 0902+343 C12 


3.39 


46.961°^} 


Q1422+231 D81 


3.10 


48.8ll°i^ 


SSA22a-MD46 


3.09 


46.761°:^? 


SSA22a-D3 


3.07 


49.7ll°f 


DSF2237+116a C2 


3.32 


/1 7 70+0.25 
^' ■ ' <J-0.25 


B2 0902+343 C6 


3.09 


45.22l}j^ 


MS1512-CB58 


2.73 


47.491}:^^ 



